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ABSTRACT 

We present the results of the first XMM-Newton observation of the interacting type Iln 
supernova 1995N, performed in July 2003. We find that the 0.2-10.0 keV unabsorbed 
flux dropped at a value of ~ 1.8 x 10""'^^ erg cm"^ s~^, almost one order of magnitude 
lower than that of a previous ASCA observation of January 1998. From all the available 
X-ray measurements, an interesting scenario emerges where the X-ray light emission 
may be produced by a two-phase (clumpy/smooth) CSM. The X-ray spectral analysis 
shows statistically significant evidence for the presence of two distinct components, 
that can be modeled with emission from optically thin, thermal plasmas at different 
temperatures. The exponent of the ejecta density distribution inferred from these 
temperatures is n ~ 6.4. From the fluxes of the two spectral components we derive an 
estimate of the mass loss rate of the supernova progenitor, M ^ 2'x 10"'' Mq yr"'^, at 
the upper end of the interval exhibited by red super-giants. Coordinated optical and 
infrared observations allow us to reconstruct the simultaneous infrared to X-ray flux 
distribution of SN 1995N. We find that, at ~ 9 years after explosion, the direct X-ray 
thermal emission due to the wind/ejecta interaction is '--^ 5 times larger than the total 
reprocessed IR/optical flux. 

Key words: supernovae: general — supernovae: individual: SN 1995N — infrared: 
stars — ultraviolet: stars - X-rays: stars — X-rays: individual: SN 1995N 



1 INTRODUCTION 



previous evolutionary stage. The interaction generates for- 
ward and reverse shock waves which bound the shocked 
wind and ejecta. This mechanism appears to be the non- 
relativistic analogue of the blast wave interaction occurring 
in 7-ray burst afterglows. The pressure and temperature be- 
hind the shocks are sufficiently high that the post-shock 
ejecta and circumstel lar material (CSM) may become p ow- 



Among all core-collapse supernovae, the particular subclass 
of Type Tin supernovae shows prominent multi-component 
hydrogen emission lines in their optical spectra. Type Iln 
supernovae have been termed also compact supernova rem- 
nants or Seyfert 1 impostors, because of the peculiar profiles 
of the emission lines and the similarity of their spectra to 
those of the Broad Line Region in AGNs. These lines are 
thought to originate from the reprocessing of radiation gen- 
erated by a violent collision of the supernova ejecta with a 
dense surrounding gas released by the progenitor star in a 



the same time synchrotron radiation is generated by elec- 
trons accelerated up to relativistic energies at the shock 
front. 




* Based on observations collected at ESO (programs 71. D- 
0265(A) and 71.D-0265(B)). Based on observations obtained with 
XMM-Newton (Obs. ID 0149620201), an ESA science mission 
with instruments and contributions directly funded by ESA Mem- 
ber States and NASA. 



f E-mail: zampieri@pd.astro.it 
I E-mail: mucciarelli@pd.astro.it 



In the "standard" model of X-ray emission from cir- 
cumstellar interaction the forward shock produces a hot 
shell (~ lO'' K), while the reverse shock produces a denser, 
cooler (~ 10^ K) shell with much higher emission mea- 
sure, from which the obser ved soft X-ray emission arises 
JChevaUer fc FranssorJll994) . The post-shock material can 
be a more efficient radiator if the explosion o ccurs in a dense 
circumstellar shell with density ~ 10^ cm~^ iTerlevich et alJ 
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Table 1. Log of the XMM-Newton and VLT observations of SN 1995N 



Start Date 


MJD 


Observatory 


Instrument 


Range/Filter/Grism 


Program/Obs. ID 


2003-07-27 


52847.2 


XMM-Newton 


EPIC pn-fMOS, OM 


0.2-10 keV, ubv UVWl UVW2 


0149620201 


2003-07-30 


52850.0 


VLT 


FORS2 


BVRI, grism 300V-f GG435 


71.D-0265(A) 


2003-07-30 


52850.1 


VLT 


ISAAC 


JHK 


71.D-0265(B) 



Il992l) . Alternatively, if the progenitor stellar wind is clumpy, 
the interaction of the forward shock with clumps of gas can 
give rise to cooler X-ray emission, with a temperature lower 
than that exp ected from interaction with a smooth wind 
JChugailllQQa) . 

Only ~ 20 supernovae have been detec ted at X-ray en- 
ergies. Recent det ection s include SN 1988Z ( Arctxaaa et^l 
1999), SN 1995N JFox et al. 200Q). SN 1998bw (Iwamoto al, 
1998), SN 1999em and SN 1999gi iSchleeel et alj|200ll) . Al- 
though soon after the explosion the X-ray emission may arise 
from Compton down-scattered 7-ray photons produced in 
the radioac tive decays of ^®Ni and ^®Co (see e.g. the case of 
SN 1987A;|ltohan|l983), in most cases evidence has been 
gathered in favour of the circumstellar interaction origin of 
the observed radiation. SN 1988Z and SN 1995N are Type 
Iln supernovae with a very large inferred X-ray luminos- 
ity {Lx ~ 10"*^ erg s"^). SN 1999em and SN 1999gi are 
"regular" Type IIP (plateau) supernovae whose X-ray emis- 
sion is about two orders of magnitude lower than that of 
Type Iln supernovae. This is usually attributed to the fact 
that in these supernovae the GSM is less dense and hence 
the ejecta interaction less strong. On the other hand, the 
powerful X-ray emission of the "hypernova" SN 1998bw is 
probably originated by synchrotron radiation produced by 
non-thermal electrons accelerated at the shock front. 

Several Type Iln supernovae show also an extraordi- 
narily large late-time infrared (IR) emission whose origin is 
still debated. The data are interpreted as thermal emission 
from dust, but its origin and the mechanism through which 
it is heated at the observed te mperatures (^ 1000 K) does 
not appear to be established fiGerardv et alj|2002ll . Several 
possibilities have been discussed, including dust forming in 
the ejecta, pre-existing dusty circumstellar medium heated 
by the supernova shock (or a precursor), IR echo powered 
by emissi on from sho ck interaction with the circumstellar 
medium (IPozzo et alJ lioOdl . For this type of supernovae 
multi- wavelength observations may provide crucial informa- 
tion about the mass-loss history in the late evolutionary 
stages of the progenitor star, the distribution of the ejecta 
and circumstellar wind (CSW), and the energy of the explo- 
sion. 

Of special interest in this context is the cas e of SN 
1995N, discovered in May 1995 JPoUas et alJll995ll . Super- 
nova 1995N is hosted in the irregular galaxy complex known 
Eis Arp 261, at a distance of 24 Mpc, whose center is occupied 
by the (IB(s)m pec) pecuHar galaxy MCG-02-38-017. The 
epoch of explosion is not known but may be e stimated to be 
about 10 months before the optical discovery llBenetti et alJ 
L1995). Franssori c t aL (^002) assumed an explosion date of 
1994 July 4, about 2 years before the first X-ray observation. 
Throughout this paper we adopt this assumption for the 
explosion epoc h. SN 1995N has been extensively observe d 
in the optical llFransson et alJl2002l : |Pastorello et alJl2005l) . 



From an analysis of the emission line profile in the optical 
and UV spectra at epochs between 321 and 1799 days af- 
ter the explosion, Franssori ct al. ( 2003) found evidence of 
three distinct velocity components for the gas, suggesting a 
clumpy circumstellar medium or an aspherical distribution 
of the surrounding gas. 

X-ray emission has been detected with ROSAT and 
ASCA at different epochs. The first X-ray observation of 
SN 1995N was perfo rmed with ROSAT HRI on 1996 July 
23 (1.3 ks;|L ewin et al. 1996..) . followed by other two expo- 
sures taken on 1996 August 12 (17 ks) and on 1997 August 
17 (19 ks; see lFox et aLlbOOOl) . An ASCA observation was 
performed on 1 998 January 19 (~ 90 ks for SIS and GIS; 
iFoxet aljlioool) . 

The high X-ray luminosity {Lx ~ 10*" - lO"^ erg s"^) 
places SN 1995N in a small group of Type Iln supernovae 
with strong circumstellar interaction (along with SN 1978K, 
SN 1986J, SN 1988Z), making it an ideal target for study- 
ing the X-ray spectral evolution and light curve of Type 
Iln supernovae. The X-ray light curve from 2 to 3.5 years 
after explosion suggests that the GSM is distributed inho- 
mogeneously and that t he average X-ra y luminosity does 
not decline significant l v l|Fox et al.ll20o61) . In the model of 
iGhevalier Sz FranssonI lll994l) a decline in flux is expected if 
the forward /reverse shock luminosity is dominated by free- 
free emission, whereas a radiative reverse shock front leads 
to a roughly constant luminosity when line emission is dom- 
inant. Both behaviours have been observed (i n SN 1986J 
and SN 1978K, respectivelvl. iFox et all ll2000l) found that 
the ASCA spectrum of SN 1995N is consistent with a sim- 
ple bremsstrahlung model with absorption {kT ~ 9 keV, 
Nh — 10^^ cm~^), although also a power-law with photon 
index a ~ 1.7 and a slightly higher column density pro- 
vides an acceptable fit. The spectral fitting hints towards 
the presence of spectral features near ~ 1.8 keV, appropri- 
ate for fluorescent silicon emission (Si K-shell lines). 

In this paper we present the results of a XMM-Newton 
observation of SN 1995N, performed in July 2003, along with 
those obtained from a simultaneous optical/infrared obser- 
vation carried out at ESO. The paper is organized as fol- 
lows. In 321 w6 present our J^MM- Newton observation and 
the adopted X-ray data reduction procedure, devoted 
to the X-ray spectral analysis of SN 1995N and two other 
field sources, while in f0]we present the OM images of the 
field. reports the X-ray light curve of SN 1995N from 
all the available X-ray measurements. In S|Elwe present the 
results from our coordinated optical/IR observations of SN 
1995N. Finally, results are discussed in ^ and conclusions 
summarized in ^ 
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Figure 1. XMM EPIC image of SN 1995N taken on July 27-28, 2003. The field is centered on the supernova and contains other two 
sources of comparable brightness, 2 and 7 (also labeled X-1 and X-3; see Table|3J. The latter is coincident with the nucleus of the galaxy 
MCG-02-38-017. [This Figure is available as a separate jpg file]. 



Table 2. Field sources (in order of decreasing RA) in the XMM EPIC image of SN 1995N 



Source number 


Source name 


RA [J2000] 


DEC [J2000] 


MOS count rate" 


pn count rate"^ 






(hh:mm:ss.sss) 


(« : ' : ") 


(10-2 s-i) 


(10-2 s-l) 


1 




14:49:08.616 


-10:18:58.75 






2 


X-1 


14:49:15.088 


-10:11:16.26 


1.15±0.11 


3.66±0.17 


3 




14:49:16.410 


-10:16:01.28 






4 




14:49:23.522 


-10:14:02.68 






5 


X-2/SN 1995N 


14:49:28.359 


-10:10:14.52 


0.96±0.07 


3.3±0.16 


6 




14:49:30.762 


-10:10:33.77 






7 


X-3 


14:49:32.720 


-10:09:48.75 


0.85±0.10 




8 




14:49:53.662 


-10:06:22.50 






9 




14:49:53.927 


-10:15:16.25 






10 




14:49:56.098 


-09:59:41.23 






11 




14:50:08.673 


-10:11:45.12 






12 




14:50:11.622 


-10:13:23.64 







" Background subtracted 



2 X-RAY OBSERVATION AND DATA 
REDUCTION 

We observed SN 1995N with XMM-Newton on July 27-28, 
2003, ~ 9 years after the estimated date of explosion (Table 
[TJl. The observation (ID 0149620201) lasted 72 ks, but was 
heavily affected by solar flares. The three XMM EPIC cam- 



eras (MOSl, MOS2 and pn) operated in Prime Full Window 
mode with the thin fllter. The EPIC MOS observation was 
split into two exposures of ~ 5 and ~ 56 ks, respectively, 
while only one uninterrupted EPIC pn exposure of ~ 64 
ks was taken. The two RGS instruments show insufRcient 
counts to perform any statistically significant analysis. 
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Table 3. Parameters of the joint fit of the XMM EPIC MOS and pn observation of SN 1995N (0.2-10.0 keV). 



Model 



JVff 
10^1 cm- 



Parameter(s) 



x2^d(dof) Flux" 



POWER-LAW 
MEKAL'' 

POWER-LAW-I-MEKAL'' 
MEKAL^'-fMEKAL*' 



n+0.4 
-0.4 

0.1+0-3 



1.4 



0.3 
-1-0.4 

0.4 



-1-0.2 
0.2 

to.9 keV 

--I-0.3 
-0.2 



1.3 



-1-0.5 
0.4 



r=2.o 

fcT=4 

r=i.5 

fcT=0.8to;l keV 



1.18(58) 
1.21(58) 
0.78(56) 



2.081-- 

^•"^-0.26 
-I 07+1. 05 
^•«'-0.70 



0.8I";| keV 0.76(56) 



:9.4 



0.1 
-1-22.6 



1.76 



+0.17 
0.38 



keV 



" Unabsorbed 0.2-10 keV flux in units of 10 erg cm ^ 
Solar abundance (frozen). 



XMM EPIC pn + MOS 



> 

S 0.01 



O 
O 

<D 



0.001 



0.0001 




Energy [keV] 



Figure 2. XMM EPIC MOS and pn spectra of SN 1995N in the 0.2-10.0 keV interval (pn above), along with the best fitting continuum 
model (solid lines) and the low/high temperature MEKAL components (dotted/dashed lines) for the three instruments. 



Data screening, region selection and event extraction 
were performed using the XMMSAS (v 5.4) software pack- 
age. Data were filtered using tlie good time intervals wlien 
the total off-source count rate above 10 keV C was < 0.5 
counts s~^ for the EPIC MOS and < 1 counts s~^ for the 
EPIC pn, leaving ~ 22 and 14 ks of useful data, respectively. 
After filtering, the first, shorter EPIC MOS exposure had a 
very small number of source photons (~ 25) and was then 
discarded. Although the standard criterion for filtering MOS 
data recommended by the XMM Science Operation Center 
(XMM-SOC-CAL-TN-0018) is C < 0.35 counts s-\ we de- 
cided to take a slightly higher threshold (C < 0.5 counts 
s~^) to have an acceptable statistics in the MOS data. Any 
residual contamination caused by solar fiares should not af- 
fect seriously our results because intense solar fiares typi- 



cally have C > 10-20 counts s"^ and EPIC MOS data 
have been analyzed jointly with the pn data (always below 
the recommended threshold). 

A filtered XMM EPIC image of the field is shown in 
Figure A list of sources with peak counts larger than the 
average counts of the adjacent background is reported in 
Table H The position of SN 1995N (RA=14h49m28s.359, 
DEC=-10°10'14".52 [J2000]) is within 0.6" from the 
accurate radio position by Van Dyk et al. (1996) 
(RA=14h49m28s.313, DEC=-10°10'13".92 [J2000]). Be- 
sides SN 1995N, other two sources appear to have compara- 
ble brightness, a field object (labeled X-1) positionally co- 
incident with a cataloged star (GSC2 S231011112849) and 
the nucleus of the galaxy MCG-02-38-017 (labeled X-3). 
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Table 4. Parameters of the joint fit of the XMM EPIC MOS and pn observation of field sources X-1 and X-3 (0.2-10.0 keV). 



5 



Source 


Model 




Parameter (s) 






X-1 
X-3 


POWER-LAW 
BREMSSTHRALUNG 
POWER-LAW 
BREMSSTHRALUNG 




r-2 1+°'^ 
i —^■'--0.1 

kT=3.5+l-l keV 


0.89(65) 
0.98(65) 
0.98(36) 
1.01(36) 



3 X-RAY SPECTRAL ANALYSIS 

3.1 XMM EPIC spectrum of SN 1995N 

Source counts were extracted from a circular region of radius 
20" centered on the source position. Background counts were 
extracted from a circular region of radius 40", on the same 
CCD. The net source count rate is reported in Table |21 A 
total of ^ 420 and ~ 460 photons were collected from the 
EPIC MOS and pn cameras, respectively. 

XMM EPIC MOS and pn spectra were binned requiring 
at least 15 counts per bin. They are shown in Figure|21 Joint 
MOS and pn spectral fits were performed in the 0.2-10.0 keV 
interval with XSPEC (v. 11.2.0). An overall normalization 
constant was included to minimize the effects of possible 
relative calibration uncertainties among the different instru- 
ments (typical fractional difference 5-10%). Despite the low 
statistics, the fit with single component models is not fully 
satisfactory (Xred — 1-2; see Table 1^. The best fit is ob- 
tained with a dual MEKAL'^ model, convolved with inter- 
stellar absorption (see Figure The best fitting column 
density of the interstellar medium is Nh ~ 1.3 x 10^^ cm~'^, 
while the temperatures of the two MEKAL components are 
kT ~ 0.8 and 9.5 keV (see Tablel^i. The improvement over a 
single component model is significant at the 4.1 a level. We 
tried also to perform spectral fits by varying the chemical 
composition of the two MEKAL components with respect to 
solar. Increasing or decreasing the chemical abundances up 
to ~ 30% in both components gives spectral fits comparable 
to those obtained with solar composition. 

Taking the best-fitting spectral model 
(MEKAL-I- MEKAL), the absorbed/unabsorbed fiux in 
the 0.2-10.0 keV band is F = 1.46 x 10"^Vl-76 x 10"" erg 
cm~^ s~^ (corresponding to a luminosity L — 9.4 x 10^® 
erg s~^, at the distance of the host galaxy d = 24 Mpc), 
about an order of magnitude lower than that of the ASCA 
observation of January 1998. 



3.2 XMM EPIC spectrum of sources X-1 and X-3 

The spectral analysis for sources X-1 and X-3 was performed 
as described in the previous section. The count rates are re- 
ported in Table |21 Source X-3 falls in a gap in the EPIC 
pn image and then only the MOS data were used. For X-1, 
joint MOS and pn spectral fits were performed. In both cases 
spectral fits were performed in the 0.2-10.0 keV interval with 
an overall normalization constant. As shown in Table|ll both 



^ The MEKAL model is the spectrum emitted by an optically 
thin, thermal plasma. 



an absorbed bremssthralung and a power-law can satisfacto- 
rily reproduce the observed spectrum, while a single black- 
body component is ruled out by the data {Xred = 2.4 and 
1.4, respectively). Assuming a power-law spectrum with the 
parameters reported in Table^I the unabsorbed fluxes in the 
0.2-10.0 keV band derived from the EPIC MOS cameras are 
F = 2 X lO"^'' erg cm~^ s~^ and _F = 3 x 10~^^ erg cm"^ s"^ 
for sources X-1 and X-3, respectively. From the parameters 
of the X-ray spectral fit and the existence of a likely coun- 
terpart identified with the GSC2 field star S231011112849, 
we tentatively conclude that source X-1 may be an X-ray 
binary in our Galaxy. The space density of such objects at 
this flux level is not known. As mentioned in §|21 object X-3 
is spatially coincident with the nucleus of the galaxy MCG- 
02-38-017. The best fltting power-law model gives parame- 
ters roughly consistent with emission from a low-luminosity 
AGN. 



4 OM IMAGING 

During the XMM-Newton observation of SN 1995N a series 
of 5000 s exposures of the field in the v, b, u, UVWl and 
UVW2 bands was also taken with the OM^. The exposures 
(3 per band) were astrometrically calibrated using standard 
IRAF tasks starting from the positions of GSC2 ESO field 
stars. The intrinsic accuracy of this procedure is 0.5". The 
calibrated images were then summed together. Nothing is 
visible at the position of SN 1995N in the u, b, v and UVW2 
bands. Using the detection limit count rates of the single 
images and the conversion factors reported in the XMM- 
SAS User's Guide (v. 3.0), we estimate the following upper 
limits for the flux of SN 1995N: _F„ < 7.7 x 10"^^ erg cm"^ 
s-\ Fb < 1.1 X 10-" erg cm^^ s'^ , F„ < 3.5 x 10"^^ erg 
cm-2 s-\ Fuvw2 < 2.3 x 10"" erg cm'^ s'^ 

A faint object is indeed detectable inside the X-ray er- 
ror box of SN 1995N in the summed UVWl image shown 
in Figure |3 The object is at the limit of detectability 
in each single frame. The position in the summed image 
is RA=14h49m28s.336, DEC=-10°10'13".77, within 0.4" 
from the radio position of SN 1995N. Within the errors 
of the astrometric calibration, we then identify this object 
with SN 1995N. The lack of catalogs of reference stars in 
the UVWl band prevent us from performing an accurate 
photometric calibration in this band. Thus, in order to es- 
timate a lower limit for the flux emitted by SN 1995N, we 

^ Note that the XMM OM u, b, v bands do not exactly coin- 
cide with the standard Johnson U, B and V band-passes. For 
a definition of the effe ctive wavelengths o f these filters see e.g. 
iRover et aP <2000l) and lMason et alJ i200ll) . 
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SN 1995N 



N 




Figure 3. OM image of the field around SN 1995N in the UVWl band (taken on July 27-28, 2003), showing that the supernova is well 
detected. 




SN iggsN 



Figure 4. VLT R-band image of the field around SN 1995N taken on July 30, 2003. 



use the detection limit count rate of a single image. Because 
SN 1995N is slightly above threshold, we can safely assume 



Fuvwi > 6 X IQ-^'^ 



erg cm 



5 X-RAY LIGHT CURVE OF SN 1995N 

Figure |K| shows the 0.2-10 keV unabsorbed fluxes derived 
from all the available X-ray observations of SN 1995N. The 
ROSAT fluxes were derived assuming a dual Raymond- 
Smith model (similar to the best fitting dual MEKAL 
model) with the same parameters reported in Table |3 For 



the A SCA fluxes, the original measurement from lFox et alJ 
l|200nh . averaged between the SIS and GIS instruments and 
extrapolated to the 0.2-10 keV interval, is F = 1.8 x 10"^^ 
erg cm~^ s~^. It should be noted, however, that in our XMM 
image two field sources (6 and 7) are sufflciently close to the 
supernova to be unresolved in the previous ASCA observa- 
tion. In particular, the flux emitted by one of these sources 
(7 or X-3) is not negligible in comparison with that of SN 
1995N. Assuming that it remained constant and subtracting 
it from that of the supernova, we flnd that the actual ASCA 
flux of SN 1995N should have been F = 1.5 x 10"^^ erg 
cm^'^ s~^. This is the value reported in Figure|21 For XMM, 
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Figure 5. 0.2-10.0 keV unabsorbed lightcurve of SN 1995N. Fluxes were de-absorbed and, when necessary, extrapolated in the 0.2-10.0 
keV interval. For the ROSATWRl data we adopt the count rates from Fox et al. (2000), assuming a dual Raymond-Smith model with the 
same parameters reported in Tablc|3] For the ASCA GIS-I-SIS data we use to the values reported by Fox et al. (2000), after subtracting 
the contribution from source X-3. The Chandra point is from Chandra et al. (2005). Finally, for the XMM EPIC data, the fluxes were 
calculated using the best fitting spectrum with the parameters reported in Table|3 The error on the XMM point is calculated propagating 
the uncertainty on the spectral parameters. 



the fluxes are the average between the different instruments 
(MOSl, MOS2, pn, respectively). The XMM value refers to 
the best fitting MEKAL-hMEKAL model {F = 1.76 x 10-^=* 
erg cm~^ s~^; see Tablej^. Finally, we report also a recen t 
measurement obtained with Chandra JChandraet alJl200£D . 
extrapolated to the 0.2-10 keV interval. 

It is worth noting in Figure |S| the large decrease in the 
X-ray flux between the ASCA and XMM observations, indi- 
cating that the supernova has started its X-ray decline. This 
conclusion is independent of the different spectral models 
used to calculate the XMM flux (see Table EJ, and IS consis- 
tent with the recent Chandra measurement. 



6 COORDINATED OPTICAL AND INFRARED 
OBSERVATIONS OF SN 1995N 

Coordinated optical and IR observations of SN 1995N were 
performed with the ESQ VLT telescopes (VLT-I-FORS2, 
VLT-I-ISAAC) on July 30, 2003 (Table[3. Two optical spec- 
tra of 1460 s each were obtained with UT4 equipped with 



F0RS2, using the grism 200V and the separation order fil- 
ter GG435. The spectral resolution is about 10 A and the 
wavelength range is 4200-9600 A. All data were reduced in 
the IRAF environment, either using IRAF tasks or specific 
procedures developed by our group. The two spectra were 
combined together and calibrated in wavelength using com- 
parison spectra of He-Ar and He-Ne lamps. The fiux cali- 
bration was performed using spectrophotometric standards. 
The estimated error on the flux is ~ 20%. The VLT spec- 
trum is shown in Figure Photometry in the BVRIJHK 
bands was performed on the VLT images, using a PSF- 
fltting technique. Photometric errors were estimated plac- 
ing artiflcal stars with the same magnitude and proflle as 
the SN, at positions close (few arcsecs) to that of the SN, 
and then computing the deviations of the measured arti- 
flcial stars magnitudes. All the photometric measurements 
are reported in Table |5] The R band VLT image is shown 
in Figure 2] 

After ~ 9 years SN 1995N is still well detected in all 
optical/IR bands. During the 1998-2003 time interval, the 
decline in the K and V band magnitudes has been ~ 3.8 and 
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SN 1995N (July 30, 2003) 




Figure 6. VLT/FORS2 spectrum of SN 1995N (3500-9600 A) taken on July 30, 2003 (grism 200V+GG435). 



~ 2.2 mag, corresponding to a flux decrement of a factor ~ 
30 and ~ 8, respectively. This is rather slow if compared 
to the typical radioactive decay decline of Type II SNe (~ 
1 mag/100 days for ^^Co), confirming that the IR/optical 
emission is powered by the interaction of the ejecta with 
the CSM (with radioactive decay of heavy isotopes playing 
only a minor role). A comparison of the decline rate of SN 
1995N in the X-ray and K bands shows that, between 1998 
and 2003 (epochs of the ASCA and XMM observations), 
the X-ray flux has decreased by a factor < 10, while the K 
band emission has dropped much more (a factor ~ 30). In 
2003 the X-ray luminosity exceeds that in the K band by 
approximately one order of magnitude (see Table [SJ. Thus, 
in 1998, Fx/Fk ~ 3. 

The optical spectrum is dominated by the H balmer 
lines, [O III] 4959-5007A, [O I] 6300-6364 A, [O II] 7320- 
7330 A (see Figure |SJ. In addition, weak emission features 
of He I are visible. The total flux of Ha decreases by a fac- 
tor 60 between the first spectrum of Fransson et al. 2002 
(1996 June 22) and our VLT spectrum (2003 July 30). In 
the same period, the flux of H/3 decreases by a factor 12 (i.e. 
the Balmer decrement reduces by a factor 5) . Similar to H/3 
is the evolution of the total flux for the He I and [O III] 
4969-5007Alines. However, over the same period the fluxes 
of the [O I] and [O II] doublets show very small changes, 
remaining roughly constant. 

The H emission lines in our VLT spectrum of July 30, 
2003, show an asymmetric profile, with two different com- 
ponents: a broader component with FWHM — 2670 km s~^ 
and a narrower one with FWHM = 1270 km s"\ As noted 



by iFransson et al.l i2002l) , the Ha profile has a significant 
evolution with time, showing a deficit of flux in the red 
wing compared to the blue wing at the same velocity. More- 
over, the broader emission component seems to peak at bluer 
wavelengths. This may be explained by the contribution of 
dust formed into the ejecta or pr e-existing in the CSM (e.g. 
iPozzo et al.ll2004l: iGerardv et aDIiooal . 



7 DISCUSSION 

SN 1995N is one of the few supernovae detected in X-rays at 
an age of ~ 9 years (see e.g. SN 1988Z). Until the last ASCA 
observation of January 1998, SN 1995N did not appear to 
show any significant decline in the X-ray flux. In fact, on the 
basis of the ROSAT and ASCA data alone, there was some 
evidence that from 2 to 3.5 years after the explosion the X- 
ray flux may hav e first dimmed by 30% and then brightened 
JFox et alj|200ol) . Our XMM observation shows that, at an 
age of ~ 9 years, almost one order of magnitude decrease in 
the X-ray flux of SN 1995N has occurred (see Figure (HI. 

Using a recent Chandra observation (March 2004), 
IChandra et alj ^2005^ find that, after subtracting the con- 
tamination of nearby sources from the ASCA measurement, 
the light curve in the soft X-ray band (0.1-2.4) appears to be 
consistent with a linear decline. However, considering also 
the hard X-ray emission, the evolution of the X-ray light 
cur ve appears more complex. According to IChueai lll993h 
and lFox et alJ ll2000l) . the early behavior of the X-ray light 
curve (up to 3.5 years after the explosion) may originate 
from a clumpy CSW, where variations of the clump proper- 
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Table 5. Optical and infrared photometry of SN 1995N from the ESO VLT observation of July 30, 2003 
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date 


MJD 


B 


V 


R 


/ J 


H 


K 


instr. 


30/07/03 


52850.0 


22.62 


21.84 


20.79 


21.00 - 






1 


30/07/03 


52850.1 








20.96 


19.75 


18.01 


2 



1 VLT+FORS2; 2 VLT+ISAAC 



Table 6. Infrared through X-ray flux of SN 1995N from our multi-wavelength observations of July 2003 



„a,b 


^UVW2 


^UVWl 


u 


pa 
^ B 


pa 


pa 


pa 


pa 


pa 


pa 


J^- '0-0.38 


< 0.23 


> 0.060 


< 0.035 




0.061-'' 


0.13-'^ 


0.059^ 


0.026-'' 


0.059-^ 





" Fluxes are in units of 10 erg cm^ s ^ 
[0.2-10] keV band 

" Efll'ective wavelength of the UVW2 filter: 2070 A iRover et a'L| 200l 
EflPective wavelength of the UVWl filter: 2905 A I Ro ver et alJboOl 
Efll'ective wavelength of the u filter: 3472 A liRovcr et alj|200d) 

^ Fractional error ~ 1% 



ties (e.g. density) with distance result in fluctuations of the 
X-ray flux. But such variations may also be consistent with 
a smooth CSW, where there is some degree of dishomogene- 
ity along the radial direction. In any case, the significant 
decline at an age ~ 9 years shows that the supernova may 
have finally entered the X-ray decline phase. Indeed, from 
all the available X-ray data, an interesting scenario emerges 
where a two-phase CSM contributes to the X-ray emission. 
Until ~ 3-4 years, the shock crosses a rather clumpy wind 
and the X-ray fiux is strongly contaminated by clump emis- 
sion, while later the CMS becomes smoother and emission 
is dominated by interaction with the diffuse gas. 

We find statistically significant evidence for two ther- 
mal (MEKAL) components from the best fit of the XMM 
EPIC spectrum. The temperature of the cooler phase (fcTi = 
0.8 keV), is consistent with the te mperature of the ther- 
mal p lasma behind the reverse shock dChevalier fc FranssonI 
119941) . The temperature of the hotter phase {kT2 = 9.5 keV) 
is similar to the temperature of the singl e-component spec - 
tral fit of the ASCA data performed bv IFox et"al] i200(i) . 
Thus this component, that dominates the high-energy tail of 
the spectrum, did not undergo significant spectral evolution, 
similarly to what observed for the high-tem perature com- 
pone nt of the X-ray spectrum of S N 1978K l|Schlegel et alJ 
I200I). As noted bv lFox et al]ll2000tl . the temperature of this 
hot gas is lower than expected for emission from the forward 
shock wave. However, if we interpret it as the temperature of 
the gas heated by the forward-shock, we can derive the expo- 
nent of the ejecta density distribution n from the expression 
T1/T2 = (3 — s)^/{n — 3)^ , where s is the expo nent of the 
CMS density distribution jFransson et aljll996ll . Assuming 
a constant and homogeneous stellar wind {s — 2), from the 
values of the temperatures Ti and T2 inferred from the X-ray 
spectrum we obtain n — 6.4, rather l ow but still consistent 
with t he range of values reported bv lChevalier fc FranssonI 
il994h (see also the recent ly reported value n ~ 5 for SN 
1978K; ISchleeel et ai]l2004^ . 

From the X-ray luminosity of the reverse shock it is 
possible to derive an estimate of the mass loss rate of the 
super-giant progenitor from which the CSW originated. As- 



suming that the post-shock gas is cooling efficiently through 
line emission, the total lum inosity of the reverse shock is 
JChevaher fc Franssonll200 J) 

Lr = 1.6 X IQ-^'cn^^Vf ergs"' , (1) 

^CSltl 1 

where M-5 is the mass loss rate in units of 10~^MQyr~^, 
Vcswi the CSW velocity in units of 10 km s~^, V4 the ejecta 
velocity close to the shock in units of lO'* km s~^, and 
Cn ~ {n ~ 3)(n — 4)/(n — 2)^. From the unabsorbed flux of 
the low temperature MEKAL component {Fx,i ~ 3 x 10^^* 
erg cm"^ s~^) we derive = 2x 10^® erg s~^ (d = 24 Mpc). 
Taking n — 6.4 we then obtain M-s/vcbwi — 0.13V4~^. As 
reported in the previous section, the simultaneous VLT spec- 
trum shows that the velocity of the broader, faster compo- 
nent of the H balmer lines is 1/ ~ 2700 km s~^. Taking this 
velocity as representative of the bulk ejecta velocity, we find 
M r^7 X lO"^M0yr-i for Vcsn, ^ 10kms-\ 

At an age of ~ 9 years, the X-ray luminosity of SN 
1995N is still dominated by a high temperature thermal 
component. If this represents emission from gas behind the 
forward shock, we can use also the luminosity of this compo- 
nent to derive an estimate of M. In fact, assuming electron- 
ion equipartiti on, the total free-free lumin osity of the for- 
ward shock is llChevalier fc Franssonlf200lfl 

where gff is the free- free Gaunt factor. Because in the 
typical physical conditions behind the forward shock elec- 
trons and ions are not in thermal equilibrium, equation ^ 
is meant to provide only a crude estimate of M-^/vcswi- 
From the unabsorbed fiux of the high temperature MEKAL 
component {Fx,h — 1-5 x 10"^'^ erg cm~^ s~^) we derive 
L/ = 10*° erg s"^ (d = 24 Mpc). At t ~ 3300 days (epoch 
of the XMM observation), we then have M-5/vcswi ~ 30 
and hence M ~ 3 x 1O~*M0 yr~^. The two estimates of the 
mass loss rate prior to the explosion of the progenitor of SN 
1995N appear to be roughly consistent. Averaging between 
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Figure 7. Infrared through X-ray spectral energy distribution of SN 1995N from our coordinated ESO- XMM observations of July 
2003. The absorbed/unabsorbed X-ray best fitting spectrum is shown as a dotted/dashed line. The optical photometry is also shown, 
superimposed on the spectrum (solid line). The arrows represent lower/upper limits for the XMM OM. 



them, we find M ~ 2 x IQ-^Mq yr~^ , at the upper end of 
the interval exhibited by red super-giants. 

Albeit the smooth wind/ejecta interaction model ap- 
pears to be qualitatively consistent with the XMM data, 
there are however some quantitative d ifficulties. The value of 
the c olumn density of the cool shell is dChevalier fc FranssonI 
l200lD 

N,^10'\n~4)^V,-' ( ^ ' cm-^ (3) 

^ 'vcsu,i * l^lOOdays J ^ ' 

For M-5/vcswi ~ 20 and V4 ~ 0.27, we obtain iVc ~ 5 x 10^^ 
cm"^, significantly larger than the value of Nh derived from 
the X-ra y spectral fit of the XMM dat a. Furthermore, in the 
model of IChevalier fc FranssonI il994) the forward shock is 
assumed to be adiabatic, while the reverse shock is radiative. 
The critical time when the cooling tim e equals the time since 
explosion f or the forward shock is dChevalier fc FranssonI 
11994 l200lD 

t..l.3xl0^(^)fi^)l/-.. (4) 

For typical parameters, ~ 10^ s and hence radiative cool- 
ing is not significant for the forward shock (while it may 
be important for years for the reverse shock). Using again 
M-s/wcsuii ~ 20 and V4 ~ 0.27, the critical time for SN 
1995N is tc ~ 1 year. Radiative cooling may still be signif- 
icant at the time of the XMM observation. Assuming that 
the forward shock emission dominates, Lf (x and the 



forward shock luminosity at the epoch of the ASCA obser- 
vation would have been 2.5 times larger, which is not suffi- 
cient to account for 0.2-10 keV flux inferred from the ASCA 
measurement, even taking into account the contamination 
from nearby sources. This may be further evidence that the 
emission up to 3.5 years was contaminated by clumps. It 
should be noted, however, that a certain spectral evolution 
at high energies between the ASCA and XMM observations 
would be expected, whereas the observations do not show 
any evidence of it. 

Because of the decrease in the X-ray flux and the 
strong contamination by solar fiares, the statistics of the 
EPIC spectrum is not sufficient to detect X-ra y spectral 
lines. Therefore, we can not test the suggestion bv lFox et alJ 
(I2OOOI) about the possible existence of a Silicon spectral fea- 
ture near 1.8 keV, although some residuals are present in 
the EPIC pn spectrum at 1.8, 2.6 and 5.4 keV. However de- 
cline in the X-ray fiux is signaling that the SN has started 
to evolve towards the remnant stage, in agreement also with 
the increase in strength of the forbidden O lines observed in 
the optical spectrum Qr'astorcllo et al. 2005). Future X-ray 
observations will then be important for detecting X-ray line 
emission. In this respect, we note th at, in a Chandr a spec - 
trum taken almost 9 months later, IChandra et alJ (|200^ 
report evidence of emission features around 1 keV that they 
identify as Ne lines. 

Thanks to our coordinated XMM-ESO VLT observing 
program we have been able to construct the simultaneous 
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broadband flux distribution of SN 1995N from the IR to 
the X-ray bands (see Table |S| and Figure |7J . Most of the 
total energy output comes from the direct X-ray thermal 
emission due to the wind/ejecta interaction, whereas, at ~ 
9 yeaxs after the explosion, reprocessed optical and near IR 
emission are less important. The K band emission of SN 
1995N is believed to be produced by reprocessing of X-ray 
radiation by d ust formed in the eje cta or pre-existing in the 
CSM (see e.g. lGerardv et alJlioCMl . Together, the fluxes in 
the R and K bands amount to « 15% of the X-ray flux. 
The R band flux is approximately 2 times larger than that 
in the neighbouring optical bands because of the signiflcant 
contribution of a few prominent lines (Ha, [OI] 6300-6364, 
[Oil] 7320-7330). In particular, the Ha line is thought to 
originate from recombination and collisional excitation of 
gas in the photo-ionized ejecta and in the cool shell bound 
by the reverse shock and the contact discontinuity. The Ha 
flux measured from the optical spectrum is 5.5 x 10~^^ erg 
cm^ s~^. This is ~ 20% of the flux of the low temperature 
component. This r esult does not appear to be in agreement 
with the findings of lChevalier fc FranssonI l)l994l that ~ 1 % 
of the reverse shock luminosity is emitted as Ha, fairly inde- 
pendent of density and other physical parameters. However, 
it should be noted that signiflcant reverse shock emission 
may be present in the far-UV, soft X-ray (below 0.2 keV) 
spectral region, thus increasing the total X-ray flux. 

The detection of SN 1995N in the UVWl band of the 
XMM OM shows that significant UV emission accompa- 
nies the X-ray luminosity, in agreement with the reprocess- 
ing of the X-ray flux that is expected to occur in these 
environments. The origin of this emission is usually as- 
cribed to Lya and other lines of highly ionized species, 
as C III-IV, NIII-V, O HI- VI, that are produced in the 
un-s hockcd ejecta, photo-ionized by the reverse shock (e.g. 
IChe valicr & Fransson 2001). 

At the epoch of our multi-wavelength observation, the 
IR emission has greatly diminished. The power emitted in 
the K band amounts to ~ 7% of the total X-ray flux. This 
late time near IR luminosity is consistent with the sim- 
ple picture of re-processing of a small fraction of the op- 
tical and X-ray radiation by a pre-existing dusty circum- 
stellar medium hea ted by the supernova shock (see e.g. 
lOerardv et alJl2002h . As the shock moves away from the 
region where the bulk of the dust is located, the emission 
quenches. 



8 CONCLUSIONS 

We observed SN 1995N with XMM-Newton at ~ 9 years 
after the estimated date of explosion. We found that the 
0.2-10.0 keV flux of SN 1995N dropped at a value of ~ 
1.8 X 10"^"^ erg cm~^ s~^, about one order of magnitude 
lower than that of a previous ASCA observation performed 
more than 6 years ago (January 1998). The decline in the 
X-ray flux is signaling that SN 1995N has probably started 
to evolve towards the remnant stage. 

Interpreting the evolution of the X-ray light curve is 
not straightforward. A complex scenario where a two-phase 
(clumpy and smooth) CSM contributes to the observed X- 
ray emission is consistent with the available data. 

The EPIC spectrum of SN 1995N shows statistically 



significant evidence for the presence of two distinct thermal 
(MEKAL) components at temperatures of 0.8 and 9.5 keV, 
respectively. If they are interpreted as the temperatures of 
the gas behind the reverse/forward shocks, we derive that 
the exponent of the ejecta density distribution is n = 6.4. 

Assuming that the observed X-ray luminosities of 
the hard/soft spectral components represent the for- 
ward/reverse shock emission, we derive an estimate of the 
mass loss rate of the progenitor, M ~ 2 x 10~^Mq yr~^, at 
the upper end of the interval exhibited by red super-giants. 

The X-ray data were obtained within the framework 
of a multi-wavelength (IR/optical/UV/X-ray) observational 
campaign. The simultaneous infrared to X-ray fiux distri- 
bution of SN 1995N shows that the direct X-ray thermal 
emission due to the wind/ejecta interaction is ~ 5 times 
larger than the total reprocessed IR/optical fiux. At the 
epoch of our multi-wavelength observation, the IR emission 
has greatly diminished. 



ACKNOWLEDGMENTS 

This work has been partially supported by the Italian Min- 
istry for Education, University and Research (MIUR) under 
grant PRIN-2002-027145. We thank an anonymous referee 
for valuable comments. 



REFERENCES 

Aretxaga, I. et al., 1999, MNRAS, 309, 343 

Benetti, S. et al., 1995, lAU Circ. 6170 

Chandra, P., Ray, A., Schlegel, E.M., Sutaria, F.K. & 

Pietsch, W. 2005, ApJ, in press 
Chevalier, R. & Fransson, C. 1994, ApJ, 420, 268 
Chevalier, R., & Fransson, C, 2001, in "Supernovae and 

Gamma-Ray Bursts", edited by K. W. Weiler (Springer- 

Verlag) ( astro-ph/01 10060 ) 
Chugai, N.N. 1993, ApJ, 414, LlOl 
Fox, D.W. et al., 2000, MNRAS, 319, 1154 
Fransson, C. et al. 1996, ApJ, 461, 993 
Fransson, C. et al. 2002, ApJ, 572, 350 
Gerardy, C.L. et al. 2002, ApJ, 575, 1007 
Itoh, M. et al., 1987, Nature, 330, 233 
Iwamoto K. et al., 1998, Nature, 395, 672 
Lewin, W.H.G. et al.,1996, lAU Circ. 6445 
Mason, K.O. et al. 2001, A&A, 365, L36 
Pastorello, A. et al. 2005, in preparation 
PoUas, C. & Albanese, D. 1995, lAU Circ. 6170 
Pozzo, M. et al. 2004, MNRAS, 352, 457 
Royer, P. et al. 2000, A&A Suppl. Ser., 145, 351 
Schlegel, E.M., 2001, ApJ, 556, L25 
Schlegel, E.M. et al. 2004, ApJ, 603, 644 
Terlevich, R. et al., 1992, MNRAS, 255, 713 
van Dyk, S.D. et al., 1996, lAU Circ. 6386 

This paper has been typeset from a TJ5X/ I^TJrjX file prepared 
by the author. 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/05 1 0046v 1 



